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Abstract 

We prove the existence and uniqueness, for wave speeds sufBciently large, of monotone trav- 
eling wave solutions connecting stable to unstable spatial equilibria for a class of A/^-dimensional 
lattice differential equations with unidirectional coupling. This class of lattice equations includes 
some cellular neural networks, monotone systems, and semi-discretizations for hyperbolic con- 
servation laws with a source term. We obtain a variational characterization of the critical wave 
speed above which monotone traveling wave solutions are guaranteed to exist. We also discuss 
non-monotone waves, and the coexistence of monotone and non-monotone waves. 

1 Introduction 

Let Z"^ be the D-dimensional integer lattice, and let p := (pi, . . . ,pd) be a point of the lattice. We 
choose and fix a vector ^ = (^i, . . . , ^/j) e such that ||^|| = 1. We define a map 

X : ^ (Z^)^+i, x(p) = (P, Xib), • ■ • , Xn{p)) 
satisfying the following: first, the translation-invariance condition 

x{p + q) = xip) + q 

holds for all p,q £ Z^; second, all points Xi{p) ^ lie on one side of the hyperplane through p 
determined by ^: 

^(p-x.b))>o, {!,..., iv}. (1) 

We shall refer to this latter condition as unidirectionality (with respect to ^). If we take D — 2, 
N = 3, and ^ = (1, 1), for example, the map 

X((fc, j)) = ((fc,j), {k - 1, j), {k,j - 1), (fc - 1, J - 1)) 

satisfies both the translation-invariance and unidirectionality conditions. 

In this paper we study differential equations on with first-order unidirectional coupling: to 
each p e Z^ we associate a function Up(-) : M M, and consider equations of the form 

Wp(i) =g(wp(i),Uxi(p)W, (2) 

Equation ([2]) is an example of a lattice differential equation (LDE). LDEs arise in diverse appli- 
cations and have been extensively studied (see, for example, ^ [51 HI [T7| [231 [T] and the references 
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therein). The unidirectional coupHng in ^ — that is, the requirement that aU of the lattice points 
determining u'p{t) lie on one side of a hyperplane through p — admittedly makes ^ into a spe- 
cialized class of problems. Such equations do, however, arise: in addition to models where spatial 
interactions are genuinely discrete and unidirectional, ^ includes discretizations of certain first- 
order PDEs and models in which the lattice sites represent agents or categories (rather than spatial 
locations). We discuss such examples at greater length below in Section [O] 

A traveling wave solution of speed c in direction ^ of ([2]) is a solution that satisfies, for all lattice 
points p and times t, 

Up{t)^q){^-p~ct), (3) 

where c > 0, ^ is as in ([T]), and the wave profile (/) is a function from M to M. We shall adhere to this 
notation throughout, and shall refer to a traveling wave solution of speed c in direction ^ simply as 
a traveling wave solution. 

All instances of ^ we consider here will have the equilibrium solutions 

Up{t) = 0\f t,p and Up{t) = l\f t,p. 

If the wave profile for a traveling wave solution satisfies the boundary conditions 

0(x) — 1 as X — — oo and 4>{x) — > as a; — > oo 

we shall say that the solution connects 1 to (though observe, in this case, that Up{—oo) = and 
Up{oo) = 1 for all p G 2.^). If 4>{x) is strictly decreasing for all a: G K, we say that the solution is 
decreasing. If (/'(x) is strictly decreasing for all x less than some Xq G K, we say that (j) is initially 
decreasing. For brevity, we shall often refer to traveling wave solutions connecting 1 to as fronts. 
Abusing terminology, we shall refer to the corresponding wave profile as a front as well. It is the 
fronts — the traveling wave solutions connecting 1 to — that occupy most of our attention in this 
paper. 

If : M — > M is the wave profile for a traveling wave solution of ^ , then (j) solves the wave profile 
equation 

C(/)'(x) = -g{(j){x), (j){x - Ki), ...,(l){x- kn)), (4) 

where Ki = ^ ■ {p — Xiip)) for i = 1, . . . , N. Conversely, if (/) satisfies Q, then Up{t) — 0(f ■ p — ct) 
satisfies ([2]). The unidirectionality condition ([T]) implies that each Kj is nonnegative, and so Q is a 
delay equation (rather than an equation of mixed type). The relative tractability of Q in this case 
is the primary reason that we consider unidirectional models. By modifying the function g and the 
constant N we can, without loss of generality, assume that > for i > 1, and we do so henceforth. 
If a particular function (/) : K — M is understood, when convenient we shall write 

^{x) = {c^{x),(t){x - Ki), . . . , 0(x - kn)) e M^+\ 

and shall write Q in the more compact form 

c0'(x) = -g(<&(x)). 

Since we are assuming that the spatial equilibria and 1 solve ([2]), the constant functions and 
1 are equilibrium solutions of Q as well. The traveling wave solution connects 1 to if and only if 
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(j) satisfies the boundary conditions 



- 1; (5) 
(j){oo) = 0. (6) 

In this paper we study the existence, uniqueness, and monotonicity of traveling wave solutions 
of ([2]) connecting 1 to 0. This amounts to considering the existence, uniqueness, and monotonicity 
of solutions of the boundary value delay equation ( 4| 5| 6 1. (Since any time-translation of a traveling 



wave solution is still a traveling wave solution, by uniqueness we shall always mean uniqueness up 
to time translation — equivalently, uniqueness of 4> up to translation in x.) Our basic approach 
for existence and uniqueness is to exploit the one-dimensionality of the unstable manifold at 1 of 
Q, denoted W^{1). We show that, for c sufficiently large, the initially decreasing branch of 
remains monotone forever, and is in fact a front. Using an Implicit Function Theorem argument, 
we show that the set U of wave speeds c for which fronts exist is open and, moreover, that some of 
the corresponding fronts are non-monotone. We also explore the behavior of waves as c approaches 
and reaches the boundary of Z^. 

To make use of results on invariant manifolds for Q, we shall need to work in the standard 
phase space for delay equations: with r — maxj(Kj), we write C — C([— r, 0],M) for the Banach 
space of continuous real-valued functions from [— r, 0] to M, equipped with the sup norm. If y is any 
continuous function defined on [r — r, t] , in the usual way we write j/r for the member of C given by 

J/t(s) = y{T + s), s £ [-r,0]. 

Given yQ G C, our hypotheses on g will guarantee that yo extends uniquely to a function y : [— r, T) — > 
M such that y solves Q for all t > and T < 00 is maximal (see, for example. Chapter 2 of ^11 J. 

We now describe the conditions that we impose on the feedback function g in ([2]). We write the 
general point of as s = (sq, si, . . . , sat), and for a € M we write a = (a, a, . . . , a) G M^"*"-^. 

The reader may find it helpful to keep the following example in mind, which furnishes a primary 
motivation for this paper: _D = l,iV = l,^ = l, and, for all p G Z, 

-Up = -Up + 2up_i - (7) 

([7]) arises in a clock model for a dilute gas described in [3T]. We shall discuss this equation further 
in section [lT] In the case of ([t]), g becomes 
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^M, g{so,si) = -so + 2si~ sf. 



The following hypotheses (Gl) and (G2) serve to generalize certain features of this motivating 
function g. Hypothesis (G3) will apply to certain examples that we consider in Section [6j 

Hypothesis 1.1. g : M is a locally Lipschitz function, on neighborhoods of and 1. g 

satisfies the following: 

• (Gl) 

- (Gl.l) 5(0) = 5(1) = 0; furthermore g(a) ^ for any a G (0, 1). 
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(G1.2) |2_(a) < for a = 0, 1; ^(a) > for i = 1, . . . , and for a = 0, 1; and 



i—O i—O 

- (G1.3) There is some /3 G with X;^o ft > 0, A > for i = 1, . . . , iV, and ft < 0, 

such that < g{s) < X^ilo ft'^* whenever 

< So < Si < 1 for alH = 1, . . . , iV. 

• (G2) If 1 < mini Sj and sq — max^ s^, then g{s) < 0. 

• (G3) There are constants Mq > and ry e (0,1) such that, whenever |so| = max; \si\ and 
\so\ > Mo, then Mfl+M < ^. 

Conditions (Gl) restrict the feedback function g in the cube [0, 1]^+^. (Gl.l) ensures that and 
1 are equihbria of ^ and Q . It further ensures that there are no other equihbria between and 1 
at which a wave emanating from 1 might get stuck. We wiU show below that (G1.2) guarantees that 
1 is a hyperbohc equihbrium with a one-dimensional unstable manifold and that, for c large enough, 
the dominant eigenvalue for the linearization of Q about is real and negative. (G1.3) guarantees 
that any initially decreasing solution of ( 4| 5 1 remains strictly decreasing at least until it crosses zero 



(see Lemma 3.1); moreover, for c large enough relative to /?, such a solution cannot overshoot (see 
Proposition 3.2). All the conditions (Gl) together will thus guarantee that, for c large enough, the 
initially decreasing branch of the unstable manifold at 1 is in fact a strictly decreasing front. The 
second inequality in condition (G1.3) can always be satisfied with some /3; the idea is, given g, to 
choose /3 to yield maximal information about the range of c for which decreasing fronts exist. 

Conditions (G2) and (G3) restrict the coupling function g away from the cube [0, 1]^. We have 
given conditions on g that are easy to verify, but our application of these conditions is embodied 
in the following two lemmas; in our main theorems, (G2) and (G3) can replaced by any conditions 
that make the conclusions of Lemmas |1.2| and |1.3[ respectively, hold. 

Lemma 1.2. Suppose that (G2) holds. Then solutions of Q with initial condition 4>q G C[— r, 0] 
satisfying 

1 < 00 (s) < '/'o(O) for all s e [-r, 0] 

are strictly increasing for t > 0. 
Proof. Since 

c0'(O) = -gim.H-i^i), H-^n)) 

and 

1 < min (j){—Ki) < max(j){—Ki) < (j){0) 

i i 

by hypothesis, (G2) implies that 0'(O) > 0. Imagine that there is some first positive time t > such 
that (/)'(t) = 0. Then, since (f) is strictly increasing on [0,t], we have 

1 < min(/)(r — Ki) < max0(T — Ki) < 4>{t), 

i i 

whence 0'(t) > — a contradiction. □ 
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As we shall see, assuming that (Gl) is satisfied, (G2) is a sufficient condition for uniqueness 
of fronts because it prevents the initially increasing branch of the unstable manifold at 1 from 
approaching zero. 

(G3) ensures that any solution of Q that satisfies the boundary condition ^ and does not 
approach ±00 remains bounded by some uniform constant M, as we now explain. 

Lemma 1.3. Suppose that (G3) holds. Then there is a positive constant M such that any solution 
of Q that satisfies ([s]) and does not approach ±00 as x ^ 00 remains bounded by M . M does not 
depend on c > 0. 

Proof. Choose Mq > 1 such that max; \si\ — \sq\ > Mq implies that \g{s) + so| < ??|so|; for some 
r] e (0, 1). Now take M > Mq. 

Suppose that supj.gR \ > M, and set 

xo = inf{ X e M : \<f>{x)\ >M} 

(and so \4>{xo)\ = M). Then \(j}{xo - < Af = IH^o)] for all i and so \g{^{xo)) + (t){xQ)\ < 
r]\(j){xo)\ < \(f>{xo)\. Thus (j){xQ) and g{^{xo)) necessarily have different signs with |(7(<i>(a;o))| > 
(1 - r,)M > (1 - ?7)Mo. 

Thus |(/)| is strictly increasing at x = xq, with |0(a;o)|' > ^(1 — f])Mo. A similar argument shows 
that |(/)(xo)| is strictly increasing, with derivative bounded away from zero, whenever xq is the first 
time that |(/'(a;o)| = M, for any M > M. We conclude that, if \4>{x)\ is not bounded by M, then 
0(00) = ±00. □ 

To complete the preparation necessary to state our main theorems, we define some quantities 
related to the linearizations of Q about its equilibria. Let c > and a e be given, with 

Ui > ioT i — 1, . . . , N . Recall that k = . . . , kn) with > for z = 1, . . . , iV is the vector of 
delays in equation Q. We define the following function of the complex number A: 

N 

D{X;c,a) -.^ cX + ao + ^aiC"'''^. (8) 

With a = {a,..., a) equal to or 1, D{X;c,Vg{a)) = is the characteristic equation for the 
linearization of Q at the equilibrium a. 
With a as above, we define the number 

c(a) inf "° + ^'f (9) 

A>0 A 

c(a) is finite if X]i!=o '^i > 0' particular, (Gl) implies that c{\7g{0)) and c(/3) are finite. As we shall 
see below, for c > c(Vg(0)) the dominant root of the characteristic equation of Q at is real and 
negative. {c{a) is the so-called "spreading speed" for the linear equation u'p — aoUp + X^iLi '^i'"xi(p) 
— see [31].) We will show that decreasing fronts necessarily exist for c > c(/3), where f3 is as in 
(G1.3). 

We also define the number 

b{a) :— inf{c > : all roots of D{\; c, a) have negative real part if c > c }. (10) 
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As c drops below 6(Vg(0)), the equilibrium solution of Q becomes unstable. 

We now state our main results. We prove our results and present examples in subsequent sections. 
We remind the reader by a "front" we mean a traveling wave solution of ([2| connecting 1 to or a 
solution of ( |4|5j6l ). 



Theorem 1.4 (Existence and uniqueness of initially decreasing fronts). Let U he the subset of 
c e [b{Vg{0)),oo) for which an initially decreasing front exists, and let Ai c U he the subset of 
c G [6(V(7(0)), oo) for which a strictly decreasing front exists. 
Assume that (Gl) holds. Then 

1. < &(Vg(0)) < c(V.g(0)) < c(p). 

2. For every c Cz U there is exactly one (up to translation) initially decreasing front. If (G2) 
holds, this front is the unique front (not just the unique initially decreasing front). 

3. Ai is closed and contains a maximal ray [cm, oo), where Cm < c(/3). Moreover, Ai does not 
intersect [h{V g{0)) , c{V [g{0))) . 

4. If g is everywhere, U is open relative to \b{y g{Q)),oo), and contains a maximal ray with 
infimum Cf, where Cf < Cm- 



Remark 1.5. In the case /3 = Vg{0), part 3 of Theorem 1.4 implies that Cm — c{j3). Part 4 implies 
that, if g is C^, there are wave speeds for which non- monotone fronts exist. We suspect that this 
smoothness assumption can be relaxed. 

Theorems asserting the existence of a critical wave speed c* such that some monostable evolution 
equation has monotone traveling wave solutions of all speeds greater than c* have a history going 
back to the 1930s [101 [T5]. We do not attempt to survey this literature here. We do mention, 
however, the distinction between the cases Cm = c(Vg(0)) and Cm > c(Vg(0)), which is discussed 
at length in [30]. In the former case, the front 0c„ travels at the spreading speed of ^ linearized 
about zero. For this reason it is called a "pulled front" because it appears to be pulled by its leading 
edge. In the latter case, 0c„ is called a "pushed front" because the speed is determined not by the 
leading edge of the front but by the behavior of feedback function g in the interfacial region. In our 
framework, we necessarily get pulled fronts when we can take /3 = \/g(0), but when we cannot take 
/3 = V(7(0), it might be that only pushed fronts are possible. 

Our second main theorem describes the behavior of fronts as c approaches and reaches a point 
of dU {cf in particular), provided that this point is greater than b{Wg(0)), that g is C^, and that 
(G3) holds. 

Theorem 1.6. Suppose that g is and that (G3) holds. If c G dU and c > &(V(7(0)), Q has 
two distinct hounded nonconstant solutions andip^, neither of which is connects 1 to 0, such that 
tP\-oo) = 1 and V/(oo) = 0. 

Intuitively, the above theorem describes the splitting, as c reaches dU from the interior of U, 
of a front — a heteroclinic between equilibria — into a pair of heteroclinics (the solutions ■0° and 
•0^) linking and 1, respectively, to some other invariant set. Under some additional hypotheses, 
we can make this idea more precise. The following, for example, is an immediate consequence of 
Mallet-Paret and Sell's Poincare-Bendixson Theorem for delay equations (see Theorem 2.1 in 
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Proposition 1.7. Suppose that N ^ 1, that g : — > M is and everywhere strictly decreasing in 
its second coordinate, and that the equilibria of Q are all isolated. Then given any bounded solution 
^ o/ Q (in particular ) either the w-limit set of (j) is a single periodic solution, or the u- and a- 
limit sets of any solution in the uj-limit set of (j) consist of equilibria of Q. 

An open problem is what conditions on g guarantee that cj — 5(Vg(0)) — that is, whether 
(roughly speaking) fronts persist as far as is suggested by the linearization of Q about 0. When 
D = 1^ N — \, and g satisfies appropriate hypotheses (such as monotonicity, and conditions to rule 



out equilibria other than 1 and 0) Theorem 1 1 .6| together with Proposition 1.7 reduce the proof that 
Cf = b{yg{0)) to showing that for c > b{Vg{0)), Q admits neither periodic solutions that oscillate 
about zero nor solutions that are homoclinic to 1: for then the solution in Theorem 1.6 cannot 
exist for c > b{'Vg{0)). Proving nonexistence of periodic solutions oscillating about for Q when 
c > 5(Vg(0)) is obverse to the problem of proving the existence of periodic solutions oscillating 
about an unstable equilibrium. This latter problem has been extensively studied for delay equations 
in the so-called negative feedback case [261 EU HSl H] • 

1.1 Related work and examples 

In this section we describe some closely related previous work and present some applications of our 
results. The part of Theorem 1 1 .4| pertaining to the existence of monotone fronts overlaps with some 
earlier results in the case that g satisfies appropriate monotonicity conditions — for example, that 

g{so, . . . , sn) is increasing in si, . . . , sat for s e [0, 1]^^"'^. (11) 

The part of Theorem |1.4| pertaining to non-monotone waves echoes some earlier results for very 
special cases of (non-smooth) g. The novelty of our work lies in the the relaxation of monotonicity 
conditions on g, and the extension of existence results for non-monotone waves to a (relatively) 
broad class of functions g. 

In [32], Weinberger considers difference equations of the form 

u{t + l,-)^Q{u{t,-)). 

Here the map Q acts on functions u which themselves map a spatial domain to the positive reals. 



The spatial domain is Z?-dimensional and either continuous or discrete. If condition (11) holds, 
time-r maps for ^ fit into the framework of [32j : in this case, the existence of decreasing fronts 
for ^ when c > c(^) is a consequence of results in [B2]- The main focus in |32j is on the so-called 
asymptotic spreading speed of initial data supported on a compact set; the uniqueness of monotone 
fronts, and the existence and behavior of non-monotone fronts, are not directly addressed. 

The existence results in [35] use a monotone iteration technique. The same is true of [2] 
(described further below), where a monotonicity condition on g is also imposed. Recently such 
techniques have been extended, in the setting of lattice integro-difference equations, to the case 
where the nonlinearity is not necessarily monotone but satisfies conditions which are similar to our 
(Gl). In particular [20] and [16] both obtain the existence of (not necessarily monotone) traveling 
waves as well as a variational characterization of the minimum wave speed guaranteeing monotone 
fronts. Although our setting and techniques differ, our results here can be regarded as complementing 
these latter works. 

We now turn to some more specific applications of 
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Cellular Neural Networks 

In [M] Hsu and Lin consider lattice equations of the form 



u'n = -Un + 2 + X! "'ifi'^ri-e) (12) 
i=o 

(two-dimensional equations, and equations where the coupling is not unidirectional, are also consid- 



ered in [H]). The chief motivation for [Tl] lies in the case that / is piecewise linear; in this case (12) 
becomes a so-called cellular neural network (CNN). Cellular neural networks were first introduced 
in [5] to model the behavior of a large array of coupled electronic components. 

At the cost of modifying / we can assume without loss of generality that z — 0. For simplicity 



we restrict attention to the case d = 1, z = 0, and ao,ai > so that (12) becomes 

= -u„ + ao/(u„) + ai/(u„_i). (13) 
In this particular case, the conditions (Gl) rewrite as follows: 

. (Gl.l) /(O) = 0, /(I) = and /(s) ^ ^ for any s e (0, 1); 

. (G1.2) < aof'is) < 1 for s e {0, 1}, and /'(I) < ^ < /'(O); 

• (G1.3) There is some {(3o,l3i) with /3o < and /3o + /^i > such that < -sq + ao/(so) + 
ai/(si) < Poso + Pisi for < so < si < 1. 

The existence problem in |14| is formulated in terms of increasing traveling waves of negative 



speed; a change of variables is required to convert to our setting. In terms of (13), the main 
hypotheses in [TJ] are 

. 1< A- 



aa+ai ' 

• o,if{s) is strictly increasing for s e (0, 1). 

• -So + ao(/o) + ai/(so) < /3oSo + PiSi for /3o = ao/'(0) - 1 and (3i = ai/'(0). 

The first and third conditions above are analogous to our (G1.2) and (G1.3). The second condition 



above is stronger; it should be thought of as analogous to (11) and allows for the application of 
monotone iteration techniques. Under these conditions, there is some c* > such that a decreasing 
front exists for all c > c* (Theorem 1.1 in [H], reformulated for our setting). In [H] and the com- 



panion paper [15, the authors also consider ( 13 ) with a piecewise linear / for which the monotonicity 
condition above fails, but which is simple enough to admit detailed analysis. In results analogous 
to ours, the authors describe condititions under which, as c drops below a critical level, monotone 
fronts give way to non- monotone fronts that "overshoot" and oscillate about their limit at oo. 

Clock model for a dilute gas 

Models of the type ([2| also arise in situations where the lattice sites represent agents or categories, 
rather than spatial locations. One example is considered in |31] . where the authors derive the 
equation 

Wp = -Wp + wl_-^, p G Z, (14) 
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where each Wp{t) is the proportion of particles in a dilute gas with a particular collision history. We 
refer to [3T] for the details of the derivation; the model ( 14 ) is also discussed in [9l [27l [30] . Under 



the change of variables u = 1 — w, (14 1 becomes the equation we introduced just before Hypothesis 

Up = g{up, Up-i) = -Up + 2wp_i - ul_^. (15) 



This is a special case of equation (13) above. 

Similar in spirit to 15] (and a chief inspiration for the present work) is |27) . in which Peletier 
and Rodriguez consider a lattice equation with a piecewise linear nonlinearity that is intended to 
imitate equation ( 14 ) . In [57] the authors describe the behaviors of traveling waves of different 
speeds: in particular, as the wave speed drops below a certain critical value the fronts change from 
being monotonically decreasing to being oscillatory. We will discuss [57] further in Section [6j 



Note that if we write (151 in the general form 

Up = -Up + f{Up^l) 



(16) 



the conditions (Gl) become 

. (Gl.l) /(O) - 0, /(I) = 1, and /(s) ^ s for s e (0, 1); 
. (G1.2) /'(O) > 1 > /'(I) > 0; 

• (G1.3) There is some /3i > 1 such that s < f{s) < /3is for s G (0, 1). 

If / is concave down on (0,1) we can 



In this case we can take /3o = ~1 and /3i > ^sup^g^g ^ 
take /?! = /'(O). We revisit ^ in Section |6] 



/(*) 
) ~ 



Semi-discrete advection-reaction equations 

Genuinely spatially discrete models with unidirectional coupling do arise; given the unidirection- 
ality assumption we impose on ([2]), though, perhaps the most natural source of applications of ^ 
is the spatial discretization of certain first-order PDEs. Gonsider, for example, the PDE 

ut = -f{u), + h{u). (17) 

The so-called first- order upwind semi-discretization with step-size e for this equation is 

u'n = -^(/(Wn) - /(Un-l)) + h{Un) = ffe ("n , ) • 

Since in this case we regard the lattice points as being separated by a distance of e, the traveling 
wave ansatz in this case is Un{t) — 4>{en — ct) and the wave profile equation becomes 

ccj^'ix) = i [f{cj^{x)) - f{<j>{x - e))] - h{cf>{x)). (18,) 

In this situation, we can write the conditions (Gl) as follows. 

• (Gl.l) h{Q) = h{l) = 0, and h{s) ^ whenever s e (0, 1); 

• (G1.2) h'{l) < < h'{0), /'(I) > 0, and i/'(0) > h'{0); 

• (G1.3) /3o < 1, /3o + /3i > 0, and < eh{so) -\- f{si) - /(sq) < e/3oSo + e/3isi whenever 
< So < si < 1. 
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Observe that, for these conditions to hold, the source term h must be monostable — h{s) > for 
s G (0, 1) with h'{0) > > thus semi-discrete conservation laws {h = 0), which are studied in 

[2 El 121], are excluded from our analysis. Observe also that if / is decreasing on any subinterval of 
[0, 1], then the first inequality in (G1.3) can be violated by choosing e sufficiently small. However, 
so long as the flux / is non-decreasing on [0, 1], and the source term h is monostable, the conditions 
can be satisfied for e > as small as we wish. 

Let us assume that (Gl.l) and (G1.2) hold and that / is smooth and nondecreasing on [0, 1]. Put 
c* = sup^g(-Q f'{s) > and /i* = sup^g^g j^-j h{s)/s. Then, given any (5 > 0, for e > sufficiently 
small the choice 



I3i 



makes (G1.3) hold. For such ^ = (/3o,/3i), 



c(/3) = inf 

A>0 



and Po — h^, — f3i 



K + je'^ - l){c, + S) 
eA 



As e — >■ 0, c(/3) -> + (5. Thus, for any c > c*, (18;) admits a unique decreasing front (f>e for e 
sufficiently small. For such a (j)^, we have < |/i(0e)|/(c — c*) everywhere. Thus the (f>e form a 
bounded equicontinuous family as e — )■ 0, and so converge uniformly on compact sets to a solution 



of the continuum limit of (18:) 



(19) 



(19) is a first order scalar ODE, and so can be analyzed easily: since h is positive on (0, 1), classical 
solutions exist that connect 1 to if and only if c > c* . In other words, the monotone fronts obtained 



in Theorem 1.4 for (18:) with c > c(/3(e)) are lattice realizations of the continuous monotone fronts 



of (19) 



Note that in the case that Cm < c*, the solutions to (18 



with c G [cm , ) converge pointwise 
almost everywhere, via Helly's Theorem, to a monotone function with countably many jump 



conditions 



discontinuities. Furthermore (/)* satisfies ( 19 ) on its intervals of continuity and satisfies the boundary 

1 and (j)^(oo) 

< c < c^, 



. In other words, 0* is an entropy solution of (19). In [25] 
it is shown that entropy solutions of (19) exist if and only if sup^g^^g j^-j 



suggesting 



that lim^^QCm{e) — sup^g^g^^ This analysis highhghts the fact that the estimate c(Vf/(0)) < 
Cm < c(/3) does not determine Cm exactly. The determination of in any particular case with 
c(Vg{0)) < c{/3) requires a careful analysis of the particular problem under study. 

In addition to the continuum limit e — !■ 0, one can look at the conservation law limit /i — > 
for fixed e. In this limit it is not hard to show that the corresponding wave profiles (f)^ converge, 
uniformly on compact sets along a subsequence, to constant functions. Thus the fronts that we study 
can be thought of as bifurcating from equilibria of the conservation law u'^ — ^(/(un) ~ f{un-i))- 

Two coupling terms 



Two more equations similar to ( 12 ) are the two-dimensional equation 



and the one-dimensional equation 
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where a and /3 are positive parameters. In each of these cases g is a function from M'^ to M — though, 
in the two-dimensional equation, if the wave direction is ^ — (1,1) then the two delayed arguments to 
(? in Q collapse to a single argument. We close this section with a remark about the wave direction ^ 
in the two-dimensional equation. With g and the coupling function x viewed as given, the set of wave 



directions ^ that satisfy ([!]) is closed and connected in S*^. Theorem 1.4 applies to each of these wave 
directions, with the quantity c(/?) now regarded as a function of ^ — or, equivalently, as a function 
of K. Suppose that /? satisfies our hypotheses and let F{x\(3,k) = ^°^-=^ — . As we shall see 
below, F has a unique positive critical point and a corresponding minimum value c(/3) = c(/3; k) > 0. 
Since F is smoothly increasing in each Ki (for positive x), so is c(/3; k). As mentioned above, in the 
/? = V g{()) case c(/3) is the asymptotic spreading speed for the linearization of ([2| about zero. In 
this light the monotonicity of c(/3) in k can be interpreted to mean that ^ transports mass more 
quickly in directions for which Xi{p) li^s further behind the hyperplane orthogonal to ^. The fact 
that c is smooth in n can be contrasted with [51 [T3] where a parameter that controls whether or not 
fronts are present is discontinuous with respect to the direction in which the front travels. 
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2 Linearization of ym and the unstable manifold at 1 



In this section we present some basic facts about the characteristic equation D{\\c,a) = 0, where 
2) is as in ([S]), and present some results regarding the unstable manifold for Q at 1. Much of this 
material is well-known. 

Proposition 2.1 (Description of the roots of the characteristic equation). Let < Q < i = 
1, . . . , N and > 0, i = 1, . . . , N be given. Take c > 0. We regard 

N 

D{X;c,a) -.^^ cX + ao + ^aie^^^^ (20) 

i=l 

as a function o/ A G C. We define 

"0 + Eili a,e'^'^ 



c{a) inf 



A>0 A 



(recall ( 10 1 j. 



The following hold. 

1. The real part of any complex root of D is less than any real root of D. 

2. If X^iLo — 0; then D has one nonnegative real root and one nonpositive real root, and all 
complex roots have negative real parts. If X^iLo '^i < 0, the two real roots are in fact positive 
and negative. 
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3- IfJ2o=i ^ ^' then c{a) is finite and < 6(a) < J2iLi on'^i < c{a). 

4- If X^o^i > ^' then D has real roots if and only if c> c{a). 

Remark 2.2. When ao = — l,ai>l, and iV = 1, a complete description of the roots of D is given 
in Theorem 6.1 in [21]. 

Proof of Lemma ^lA\ We begin with the case that = for all i e {1, . . . , N}; in this case we need 
to prove the first two points of the proposition (since X]o^=i ^ ^ implies that > for some 
positive i). In this case, D[X) = cA + ao has one root, this root is positive, and the proposition 
holds. We henceforth assume that ai > for at least one i G {1, . . . , N}. 
Consider the restriction 

N 

D{x; c, a) = cx + aQ + aie"'*''^ 
1=1 

of D to the real axis. Let D'{x; c, a) = c—J2iLi OiiKie~'^^^ denote the derivative of D with respect to 
X. Note that D"{x;c,a) = J2tLi ainfe''^'^ > 0, that D{-oo) = D(oo) = oo, and that D'(-oo) = 
— CJO and D'{oo) = c > 0. Thus D decreases to a single global minimum D = D{c,a) at some 
X = x{c, a), and thereafter increases. 

We now prove part (1) of the proposition. Assume that Xq — xq + iyo is a complex root of D 
(with 2/0 7^ 0)- Writing the real and imaginary parts of the equation D ~ separately yields 

ao + cxq + Y^iLi aje""'^" cos(Kiyo) = 0; 

C2/0 + Eili a,e-'''=^o sin(Kiyo) = 0. 
If the above equations hold, we see that 

N 



D{xo; c, a) = ao + cxq + ^ a^e cos(O) > 0, 



1=1 

with strict inequality unless cos(Ki?/o) = 1 for all i e {1, . . . , N}. The second equation above rewrites 
as 

Thus at Xq we have 

^Dixo;c, a) = c - y a,K,e-^''^° = - V a,K,e-^'-« ^"H^ + 1 < 0. 
dx ^ V 'i^Vf) J 

Thus, viewing Z? as a function of the real variable x, we see that D is nonnegative and decreasing 
at Xq. Thus Xq is less than or equal to any real root of D. If xq is in fact a real root of D, then 
cos(Kij/o) = 1 for s-U * G {li • • ■ ! ^-nd the equation D = must hold both at Xq + iyo and at Xq. 
The imaginary part of D = shows that this is impossible; we conclude that xq lies strictly to the 
left of any real roots of D. This completes the proof of part (1) of the proposition. 
Since 

AT 



D{Q;c,a) = ^ a^. 



i=0 
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if X]t=o cti < we see that D has a nonnegative real root and a (not necessarily distinct) nonpositive 
real root. If X]t=o cti < 0, D has a positive root and a negative root. Part (1) now implies that all 
nonreal roots necessarily have negative real parts. We have proven part (2) of the proposition. 
We now assume that J2iLo '^i ^ ^- With this assumption, let us define the function 



F{x;a) = F{x) 



N 

1=1 



(21) 



for a; > 0. Note that F{{)+) = F{oo) ~ oo and that F is continuous; thus F is minimized at some 
positive point and its minimum value is c(a). 
We compute 



F{x) 



■ N 

.1=1 



ao 



> 



1 



N 



N 



"0 



N 



Taking the minimum over x > of the left hand side yields c(a) > X]i!=i ^i^^i- To prove part (3) 
of the proposition, it remains to show that b{a) < X^iLi Q^i'^ji that is, for c > X)i!=i '^i'^ii roots 
have negative real part. 

Let us assume accordingly that c > X]i!=i Q^i^i- this case a simple computation shows that 
D{0) > and D'{0) > 0, so that any real roots of D are strictly negative. Thus, if D has real roots, 
the complex roots must have negative real parts by part (1) of the lemma. Even if no real roots are 
present any complex root must have negative real part: for given a root xq + iyo with yo ^ 0, we 
know from above that D'(xo; c, a) < 0, whence xq is less than 0. 



It remains to prove part (4) of the proposition: still assuming that ^ 



N 



> 0, we wish to 



show that D has real roots precisely when c > c{a). Let us write D for the minimum of D over the 
reals and x for the point at which this minimum is attained (our earlier observations imply that x 
is negative). We have 



N 



which we substitute into ( 20 ) to obtain 



i=l 



N 



D{x) = ao + 'y^ai{l + XKi)e 



(22) 



(23) 



Note that 



i=l 



N 



b'{i) = -i^a.K^e-"^^^ > 0. 



Since x is defined implicitly in terms of c, we may regard (23) as defining D as a function of c 



(holding a fixed). Since x is strictly decreasing in c and D is strictly increasing in x for i < 0, it 
follows that D is strictly decreasing in c. As c approaches X^iLi ^i^i from above, x approaches 
from below and D approaches X^iLo Q^i > from below. On the other hand, as c oo, x approaches 
— oo and D approaches —oo. Therefore there is some unique c* where D — viewed as a function of 
c — crosses zero downward as c crosses c* upward. When D < 0, D has real roots. 
Our claim is that c* = c(a). To establish this we compute 

N 

-ctQ + ^ ai{K.iX — l)e'^* 



F\x) = 



1 
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and observe that 

Since c* is the unique positive root of Z?, it foUows that — i(c*) is the unique positive critical point 
of -F — that is, the global minimum of F . Thus we have 

c(a) = F{~S:{c^)). 

To see in turn that = a;(c*)), we compute 

= b{c,) ^ 



N 



AT 



-a;(c*) 



F(-i(c*)) = c*, 



where the last line follows from (22 1 and from the definition of F. □ 



The above proposition establishes all of point 1 of Theorem [L4| except that c(/3) > c(V5(0)). 
Our hypotheses on /3 imply that Pi > dg{0)/dsi for all i G {l,...,iV}, and that J2^=ol^i — 
J2^=o99{0)/ds^. It follows that F{x;f3) > F{x;Wg{0)) for all a; > (see ([28|); the desired in- 
equality follows in turn. 

Our hypotheses on g and the above proposition yield that, for every c > 0, the characteristic 
equation for the linearization of Q at 1 has exactly one positive root A, with all other roots having 
negative real parts. Accordingly, equation Q has a one-dimensional unstable manifold at 1. We 
now review some facts about one-dimensional unstable manifolds at a hyberbolic equilibrium point; 
we refer the reader to chapters 7 and 10 of llj. Recall that we are writing C — C[—r,0] for the 
phase space for Q. Following [TI] let us write 

VK"(1) = { e C : y{t) solves Q for alH < and y(-oo) = 1 }. 

Given a neighborhood V of the constant function 1 in C, we write 

V)^{yQ£ W^{1) : yteV for aU t < }. 

The following proposition is just Theorem 1.1 in chapter 10 of 11^ (along with some facts that 
emerge in the proof), adapted to our situation. 

Proposition 2.3. Let c > and let A denote the unique positive real root of the characteristic 
equation for Q at the equilibrium 1 . There is a neighborhood V about 1 in C and a linear projection 

n -.C := {ke^-, keR} dC 

such that the following holds. There is some neighborhood V of inU and a function h : V ^ C 
such that W^"(l, V^) = 1 -|- h{V), h{0) = 0, and such that h satisfies the estimate 

\\hi^)-hm > u-^\\/2. 
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We apply these facts in the following lemma. 

Lemma 2.4. Given any c > 0, up to time translation there are exactly tiuo nontrivial solutions 0+ 
and (/)_ of Q which satisfy the boundary condition ([5| . There is some Xq such that (/)_ {x) — 1 is 
strictly negative and (f>+{x) — 1 is strictly positive for all x < Xg. 

Proof. Let V be as in the above proposition. Any solution of Q satisfying ^ certainly satisfies 
(j)x & W'^{1, V) for all sufficiently negative x. For such x, let us write 4>x — ^ = Vx = h{ke^'), where 
h is as in Proposition |2.3| 
Since h is C^, we have 

yx=lke^- + f{ke^-), 

where ||/(fce'*' )|| < e||fce'^ || for any given e, provided that k is close enough to zero. The expansivity 
condition in Proposition |2.3| implies that 7 is nonzero. 

Observe that the point ke^' G C has maximum absolute value |fc| and minimum absolute value 
|fc|e~'^''. Therefore, by choosing k close enough to we have that ||/(fce^^)|| < |7fc|e~^'", whence we 
can conclude that yx is all of one sign and of the same sign as ^ke^' . Thus, for any € W^{1, V) 
sufficiently close to 1, — 1 is all of one sign, and both signs are possible. 

Thus (j)- and 0+ as described in the lemma exist. The uniqueness of 0_ and c/)-^ up to time 
translation can be proven by noting that ft, is a diffeomorphism near zero; we omit the details. □ 

Proposition 2.5. Suppose that Ck Cq ^ as k ^ 00, and that a globally defined solution 
(/)fe : M — > M 0/ Q of speed Ck exists for all k. Suppose moreover that there is some M such that 
\4'k{^)\ l£ for all k £ N and all x Cz R. Then (j}^ converges, uniformly on compact subsets o/M, 
to a solution (f)Q of Q of speed Cq . 

We emphasize that, even if the ipk all connect 1 to 0, 4>o need not do so. 
Proof. Since the bounds |(/)fe(a::)| < M and 

< Bnp{ IffWI, \s\<M} 
miUk Ck 

hold uniformly for a; G M and /c G N, the sequence <^fc is uniformly equicontinuous and pointwise 
bounded and hence, by Ascoli's Theorem and a diagonalization argument, has a subsequence that 
converges on compact sets to some limit ^q- Since the (j)^ satisfy Q, 0^ is equicontinuous and 
pointwise bounded as well and hence also converges uniformly on compact sets to some limit. This 
limit is equal to —g{^Q{x))/c for all x, but is also equal to (pQ by a standard theorem (see, for 
example. Theorem 7.17 of [28 ). It follows that 0o is a solution of Q. □ 

We close this section by establishing an estimate on the basin of attraction for ^ about that 
is locally uniform in c. Suppose that Tc{t) : C C is the solution operator for Q (with c the wave 
speed) and that Lc{t) : C ^ C is the solution operator for the linearization of Q at zero: 

cy\t) = -Vg{0)-Yit), Y{t) = {y{t),y{t ~ . . . ,y{t - kn))- (24) 

As proved above, as long as c > b{Vg{Q)) there is some Ac > such that every root of the charac- 
teristic equation at zero has real part less than —Ac. For every such c and Ac, it is standard that 
there is a Kc > such that 

\\Lc{t)yo\\ < K,e'^^'\\yo\\. 
We will need the following elementary lemma. 
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Lemma 2.6. Given ci > b{T/g{0)), e > 0, and r > 0, there is a S > such that \c — ci\ < S implies 
that 

\\L,{t)yo - L,,{t)yo\\ < e\\yo\\ 

for all t e [0, r]. 

Proof. With notation as in the statement of the lemma, write c = Ci +r]. Given j/q G C', let us write 
y for the continuation of j/o a solution of ( 24 1 with wave speed Ci and w for the continuation of 
yo as a solution of ( 24 ) with wave speed c — so 

yt = Lci{t)yo and Wt = Lc{t)yo. 

Write Ayt = -V g{'Q) ■ Y{t); we have y'{t) = Aj/f/ci and w'{t) Awt/c. For alH > we have the 
bounds 

M M 
|y'WI<-||2/*||and \nj'{t)\<-\\wtl 

where M > |Vg(0)|. 

Choose c?i > such that \\yt\\ < Ill/oil for all t £ [0, t]. 

Now, observe that \\wr — VtW < /J" ll^^'l^) ~ "^^i t^^t for all t G [0, r] we have 



\Wit)~y\t)\\ = 



< 



< 







c 


Cl 


^(wt - 


2/t) 


c 




M 

— \\wt - 
c 




M 

— W^t - 

c 





Comparing — j/f || to the solution of the ODE 

, M diMri 
u = — u-\ -. r 

C Ci(Ci+77) 

yields the desired result for \r]\ sufficiently small. 



Ayt Ayt 

C Cl 

My ^ 
ci(ci + 

diMrj I 
ci(ci + 77) 



boll, «(0) =0 



ytW 

Vol 



□ 



We now make use of two standard facts. First, Lc{t) is the derivative of Tc(i) with respect to its 
functional coordinate at 0. Second, T'c(t) is uniformly differentiable at with respect to t and c, so 
we have the following: given any r/ > 0, tq > 0, ci > b{Vg{Qi)), and 5 G (0, ci — 6(Vg(0))), there is 
some e > such that ||j/o|| < e implies that 

||Tc(t)yo - ^cWyoll < ?7l|yoll for all t e [0,to] and c e [ci - 5,Ci+S]. 

(See, for example, sections VII. 5 and VII. 6 of 8 .) The following proposition can be obtained from 
Theorem VII. 1.3 in [7] together with the fact that the sun-star calculus (see [S]) can be used to 
write Q as an ODE in a Banach space. For completeness, we include the proof here. The ideas are 
similar to those in section VII. 5 of [S]. 

Proposition 2.7. Let ci > 6(V.g(0)) be given. Then there is a S small enough that the following 
holds. There is an e > such that whenever |c — ci| < S, then the basin of attraction for Q about 
includes the ball of radius e centered at the origin. 
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Proof. As above, choose K > and A > such that 

\\L,,{t)yo\\ < Ke-''\\yo\\ 

for all t > 0. 

Choose r large enough that A'e^^"^ < 1/8. Now choose S, as in Lemma 2.6 small enough so that 

\\L,,it)yo - L,{t)yo\\ < II2/0II/8 for alH G [0, r] (25) 

whenever \c — ci\ < S. Now choose e > small enough that ||yo|| < e implies that, for all t E [0, r] 
and any c G [ci — (5, ci + (5] , 

||T,(i)M-ic(%o||<^. (26) 
At t ~ T, then, for ||yo|l < e we have 

\\Tc{r){yo)\\ < \\Tc{T){yo) - Lc,(r)(yo)|| + 1 1 ^c, (r) (2/0) 1 1 

< rc(T)(yo) - ic(T)2/o|| + l|ic(T)yo - L,^{r)yo\\ + Ke'^'\\yo\\ 

< \\yo\\/2. 

Lett > be given and write t = nr+ti where ti € [0,t). Write Tc(i)yo = {Tdh) - Ldh) + £c(ti)) Tc(t)"(?/o)- 
Since ||Tc(r)"?/o|| < 2""||?/o|| < e the estimates ([25]) and ([26| hold and we have 

||r,(%„||<(^+A-)2-"||yo||, 



which goes to zero as i — cx). 



□ 



3 Existence and uniqueness of monotone fronts when c > c{(3) 

In this section we establish the existence and uniqueness of monotone fronts when c > c(/3). In 
particular, we prove statements 2 and 3 of Theorem |1.4| Already from last section, we see that the 
only candidates for a front (i.e. a traveling wave solution connecting 1 and 0) are the solutions 4>+ 
and (f>- described in Lemma |2.4[ of these, only (/)_ might be an initially decreasing front. Indeed, 
henceforth we will be preoccupied almost entirely with (/)_ and whether it satisfies the boundary 
condition We begin with the c > c(/3) case. Our proof that (f>- is actually a decreasing front 
when c > c(/3) proceeds in two steps: we show that (f>- is strictly decreasing until it crosses zero; 
and we show that never does cross zero. We assume that (Gl) holds throughout. 

Lemma 3.1. Let ip be a solution of Q that satisfies the boundary condition and suppose that 
Xq E R is such that 4>{x) € (0, 1) for all x < Xq. Then 4>' {x) < for all x < xq- In particular, if 
4>{x) £ (0, 1) for all X, then <j} is strictly decreasing everywhere. 

Proof. It is sufficient to prove that ((>' {xq) < 0. Imagine, to the contrary, that 4''{xo) > 0. Write d 
for the minimum nonzero delay Ki. We claim that there is some xi < xq — d with (p'^xi) > and 
(f'i^i) ^ 4'{xo)- Induction then contradicts the assumption that 0(— 00) = 1. 

We prove the claim. To say that ^'(xo) > is to say that g(<i>(xo)) < 0. Since (j){xQ — Hi) G (0, 1) 
for each i it follows from the first inequality in (G1.3) that there is some i with (j){xo — Ki) < (I){xq). If 
(j)'{xQ—Ki) > 0, we're done: just take xi = XQ — Ki. li(f)'{xo — Ki) < 0, then there must be a minimum y 
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on [xq — Ki, xq) such that (t>'{y) = 0. The same argument we just made shows that (f)(y — kj) < (j){y) 
for some j; since (j) is strictly decreasing on (xq — Hi,y) we see that y — Kj < xq — Hi, and that 
(j){y — Kj) < (f>{xo — Ki). By the Mean Value Theorem, (t>'{xi) > for some xi G (y — Kj,xo — Ki). 
This proves the claim and completes the proof. □ 

Proposition 3.2. If c> c{f3), the unique initially decreasing solution o/ Q satisfying the bound- 
ary condition ^ is strictly decreasing for all time and satisfies the limit ^ . In particular, it 
monotone front. 



IS a 



Proof. Let 



/3o + Eti Ae"'" 

a;>0 X 



A — arg mm 



(The minimum exists because '^^^q Pi > — recall Proposition 2.1 ) Since c > c(/3), for some 
S G (0, A) we have 

A-5 

Consider the supremum 

Xq :— sup{x I (f>'{y) > ~A(j>{y) for all y < x}. 

Since (f){—oo) = 1 while oo) = 0, the set over which we are taking the supremum is not empty. 
We will show that xq — oo, precluding that (j) is ever equal to 0. 

By way of contradiction, imagine that xq < oo. Since ip is C^, we must in fact have 0'(a;o) = 
—A(p{xo). Comparing (p to the ODE 

tP' = -A^j, ip{xo) = (t){xo) 
ixo - Ki) < 0(xo)e'^'"' 
-3(0(2^0), (Pi^o - Ki), (f){xo - kn)) 



yields that, for all Ki 
and therefore that 



0'(xo) 



Mi^o) - Eill PiH^O - Ki) 



> 



> 



/3o + Etl 



(pixo) 



= -{A- S)(f){xo) > -A(f>{xo), 

a contradiction. Note that in the second line we have used the second inequality in (G1.3); this is 
the only place in this work where this inequality is used. 

Thus the solution is decreasing for all x and bounded below by zero. It follows that the limits 
0(00) and 4>'{oo) exist, that (j)'{oo) — 0, and that 4>{oo) G [0,1). Thus (pioo) — g($(oo)), and it 
follows from (Gl.l) that 4>{oo) ~ 0. This completes the proof. □ 
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We have established the existence and uniqueness of a monotone front in the c > c(/3) case. To 
show uniqueness among all fronts, as in statement 2 of Theorem |1.4[ we must rule out the possibility 
that ^+(00) = 0. This is the purpose of hypothesis (G2). Since 0+ satisfies ^ but is not identically 
one, there must be some G M such that (/)+(a;o) > 1 and 0+(a;o) > (p+i^o + for all s G [— r, 0]. 
In this case, though, hypothesis (G2) yields that (j>'^{xa) > 0. The conditions just articulated will 
be preserved as x moves forward; thus in fact (j)+{x) is strictly increasing for all x > xq- 

Thus we have 

Proposition 3.3. If c > c(/3), is the unique monotone front connecting 1 to 0. // in addition 
(G2) holds, is the unique front connecting 1 to 0. 

Let us write 

Cm — inf{c > : every initially decreasing front for c > c is monotone } 
(this is the c„i of the theorem.) 



We now show that M is closed and contains [c(/?), 00) as in statement 3 of Theorem 1.4 



Proposition 3.4. M is closed and contains [c(/3), 00). In particular, when c — c„j, the solution 0_ 



described in Proposition 2.4 is a monotone front. 



Proof. Let be a sequence of points in Ai with limit cq and let (p^ denote the corresponding 



monotone fronts. Let us translate each 0^ so that 4>kiO) = 1/2. We know from Proposition 2.5 that 
(pk 00 J where the convergence is uniform on compact subsets and (pQ is a solution of Q. Since 
each (j)k is monotone, 0o is monotone as well and hence has limits at a: = ±cxi; additionally the 
limits (/)o(±oo) exist and are zero. Thus 0o(±oo) = g($o(±oo)). The pinning (/>fc(0) ~ 1/2 yields 
that (j>o{—oo) € [1/2,1] and (f>Q{oo) e [0,1/2]. It now follows from (Gl.l) that (/)o(— 00) = 1 and 
(j)a{oo) = 0. Since 0o(~oo) = 1 and 4>o{x) is never greater than one, it follows that 0o = 4*- and 
that (up to translation) is the unique decreasing front. □ 

We now complete the proof of statement 3 in Theorem |1.4| by showing that M does not intersect 
[6(V<7(0)),c(V5(0))). 

Proposition 3.5. Let (f)_ denote the initially decreasing branch ofW"{l), and assume that 4>- is 
a decreasing front. Then c > c(V.g(0)). In particular, Cm > c(V(7(0)). 

The ideas in the following proof are similar to those in sections 7 and 10 of [H]. We rely on 



the fact that solutions of the linearized equation (24 1 have expansions in eigenfunctions (see, for 
example, section V.3 of [8 , and [H]). 

Proof. Let ynix) :— ■ We shall write 

^-{x) = {(f)-ix),(f>^{x - Ki), . . .,<p^{x - kn)); 

Yn{x) = {yn{x),Vn{x - Kl), . . ■ ,yn{x - Kn))- 

Observe that each function ?/„ satisfies the delay equation 

1 



cp-in) 



g{cP^{n)Y^{x)) 



^-[~Vg{0)-Y^{x) + 5^{x)] 
c 
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where Sn{x) = Vg{0) ■ Yn{x) - ^-^g{(p^{n)Yn{x)). 

Since 0_ is hypothesized to be a decreasing front, we have yn(x) E [0, 1] for all x > 0. Note also 
that, for n sufficiently large, Unix) is approximated by 

-V.g(0) •r„(2; + n) 



and so the ?/„ are an equicontinuous family for a; > 0. Therefore some subsequence of {i/n} converges 
uniformly on compact subsets of [0, oo) to some yQ. The remainder term (5„ satisfies Sn = 0{(j)-{n)) 
on the half- line x G [0,oo). Thus j/q satisfies the linear constant-coefficient delay equation 

cy'o{x)^-VgiO)-Yoix) x>0, (27) 



which is simply equation ( 24 1 restricted to a half- line 



The solution yo has an expansion in eigenfunctions. Looking at the real and imaginary parts 



of the equation D{X; c, Vg(0)) = (recall the proof of Proposition 2.1 ) shows that, if xq + iyo is a 
nonreal root of D, the only roots of D with real part xq are xq + iyo and x — iyo. Moreover, since D 
is entire its roots are isolated. Thus there are at most two dominant terms in the expansion of j/o! 
these are either oscillatory (if the corresponding roots are nonreal) or monotone (if the corresponding 
root is real). 

Since each y^ is decreasing, yo is non-increasing. There are now two possibilities. Either yo is 
identically zero, or the root corresponding to the leading term in the expansion of yo is negative and 
real. But yn{0) = 1 for all n, so yo{0) = 1 and yo cannot be identically zero. This completes the 
proof. □ 

We have proven parts 2 and 3 of Theorem |1.4| 

4 U is open 



In this section we complete the proof of Theorem 1.4 by proving statement 4. In particular, we 
show that if c G W and c > b{\/g{0)), then U contains an open interval in R about c. We use a 
continuation argument. We assume that (Gl) holds throughout, and in addition assume that g is 
C\ 

We write W^'°° = W^'^{M.) for the Banach space of Lipschitz functions on M, equipped with the 
norm 

||z|| = esssup][j|z| -I- esssupalz'l, 
and = L°° (M) for the Banach space of essentially bounded functions equipped with the norm 

||z|| = csssupkI^I. 

Given G K, we write W^'°° C W^'°° for the codimension-one affine subspace of functions z 
satisfying z(0) — v. Given z, as we done have earlier we shall write 

Z{x) = {z{x), z{x — Ki), . . . , z{x — Kn)) . 

We define the following function F : W^'°^{R) x M+ ^ 

F{z, c){x) cz'{x) + g{Z{x)). (28) 
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Note that zeros of F are solutions of Q with speed c that are defined on all of M. 

Suppose that v e (0, 1) and that (j) £ W^'°° is an initially decreasing front of speed cq. Below 
we linearize F about ((/>, cq) and show that the resulting operator is an isomorphism between W^q 
and L°°. Wq'°° is the tangent space to therefore by the Implicit Function Theorem F{-,c) 

has a zero 0c in W^'°° for each c sufficiently close to cq. This 0c is a solution of Q with speed c 
that is defined for all a; € M, is translated so that 0c (0) = i^, and is close to in the W^'°° norm. 



Proposition 4.3 establishes that 0c satisfies the boundary conditions 0c(— oo) = 1 and 0c(oo) = 
— and so 0c is a front. Thus we see that fronts connecting 1 to (no longer necessarily monotone) 
exist for an open subset of wave speeds — in particular, for some wave speeds c < Cm (and so 
non- monotone fronts connecting 1 to must exist). 

The derivative of F with respect to its functional coordinate is given by 

DiF{z, c)z = cz'{x) + \/g{Z{x)) ■ Z{x). 

Observe that the map (z,c) i— > DiF{z,c) E C{W^'°° , L°°) (where the latter space is endowed with 
the usual operator norm) is continuous. (This is where we use the fact that g is C^, rather than 
only near and 1.) Therefore, if we establish that £'ii^(0, cq) is an isomorphism from T4^q^'°° to 
L°°, the Implicit Function Theorem applies. 

Proposition 4.1. Let cq > b := b{Vg{0)) and suppose that there is some initially decreasing S 
W^'°° (not necessarily monotone) that satisfies F(0, cq) = and the boundary conditions 0(— oo) = 1 
and 0(oo) = 0. Suppose further that v has been chosen so that, whenever x < 0, 4>{x) G {v, \), 
§^^mx))<0, J2l^§M.m^))<0, and 

Vg($(a;)) • (so, • • • ,sn) < when < Si < sa for all 1 < i < N. 



Then L :~ DiF((j>, Cq) is an isomorphism from Wq'°" to L°° . 

Remark 4.2. In light of (G1.2) and the fact that g is C^, the conditions on V(7($(a;)) for x < 
can be imposed without loss of generality by choosing v appropriately. 

Proof. Since cq > 6 the characteristic equation of Q at has no roots in the closed right half plane. 
Since Cq > 0, the characteristic equation at 1 has one root in the right half-plane. Thus, it follows 
from Theorem A in [22] that the Fredholm index of L is 1 when L is regarded as an operator from 
W^ °° to Since Wq^'°° is codimension one in 1/1^^'°°, the Fredholm index of L is zero when 

regarded as an operator from VFq'°° to L°° . Thus to show that L is an isomorphism it suffices to 
show that its kernel in Wq^'°° is {0}. Suppose that w G W^'°° \ {0} and that Lw — 0. We will show 
that w;^Wo'°°. 

The expression Lw — rewrites as 

cw'{x) ^ ~Vg(<^{x)) -Wix), W{x) = {w{x),w(x- Ki),. . .,w{x - kn))- (29) 

We claim that w satisfies the boundary condition w(— oo) = 0. In fact this claim follows from the 
theory of exponential dichotomies, but for completeness we provide an elementary proof here. There 
are two cases to consider. Either w is initially monotone, or it isn't. If w is initially monotone, then 
the limits ii;(— oo) and u;'(— oo) exist and u''(— oo) — 0. Thus = Vg(l) • {w{—oo),---w{—oo)). 
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Since the sum of the entries of V(7(l) is negative, we have w(— oo) = as desired. We now consider 
the case where w is not initially monotone. There is a sequence Xk — > — oo with w'{xk) = 0. 
Choose a subsequence so that w{xk) and 'w{xk — Ki) simultaneously converge, for all z, to limits 
If* = Iimsup2._^_g2 w{x) and Wi. Then, letting fc — )■ oo in the expression = Vg($(xfc)) • W{xk), we 
obtain 

= Vg(l) • (w*, wi, . . . ,wn)- 



Imagine that w^, > 0. Then 
for alH e {1, . . . , N}, whence 



dg{l) dgil) 

-T, Wi < — -I 



i=0 



a contradiction. The case < is ruled out similarly. Thus the limit exists and the claim is 
proven. 

Since limj._>._oo w{x) = 0, the solution w lies on the unstable manifold at for the linear equation 



(29 1. Thus, from an argument similar to the proof of Lemma 2.4 except with (V(/(<i>(a;)), •) replacing 
g{-), we know that w is initially strictly of one sign. Without loss of generality we take w{x) > for 
all sufficiently negative x. Choose xq < such that w{x) > for all x < xq and also such that 

w{xo) > w{xo — r), re [0, r] 

(such an xq exists because w{—oo) — 0). Our hypotheses on \/g{^{x)) for a; < now yield that 
w'{xo) > 0. 
Now let 

:= supja; e M : | w'{y) > for all y £ [xo,x)}. 

If x^, is infinite, then w ^ Wg^'°° and we are done. Suppose that x^, is finite. Then w'(x^) = so 
Vg{^{x^)) ■ W{x^) — 0. Since w{x) is strictly increasing on [a;o, a:*], by our hypothesis w'{x) > on 
[xo,x^^,) n (—00,0]. It follows that x» > 0. Thus in particular w{x) > for a; < 0, and w ^ Wg'°° . 
This completes the proof. □ 

We now complete the proof that U is open. Suppose that cq £ U, with a corresponding front 0. 
Apply the Implicit Function Theorem with c as a parameter: for each e > there is a, S = 6{e) > 
such that whenever |c — co| < S then there is a (j)c which solves Q and which satisfies ||0c — 0o|| < £■ 
Moreover, since 0c is continuous in c we have ||0c ~ '/'oil as (5 — 0. 

It remains to check that (f)c connects 1 to if (5 is small enough. Given any e, by choosing c close 
enough to Cq we can arrange for ||(/)c — < e/2 in the W^'°° norm, and so also in the sup norm. By 
shifting (j>c and (j) (identically) so that |0(a;) — 1| < e/2 for a; < 0, we guarantee that \4>c{x) — 1| < e 
for all a; < 0. Similarly, we can shift (/>c so that \4>{x)\ < e for all x > 0. That this is enough to 
guarantee (pd^oo) = 1 and (j)cioo) = follows from the following proposition. 

Proposition 4.3. For each cq > b{Vg{0)), there is an e > such that the following hold. 

1) Suppose that |c— co| < e and thatipQ G C has a global backward continuation for Q, with wave 
speed c, such that Wipr ^ 1|| < e for r < 0. Then the limit limt_>._oo ''Pit) exists and is equal to 
1. 
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2) Suppose that |c — co| < e and that ^po G C has a global forward continuation for Q, with wave 
speed c, such that \\iPt\\ < e for r > 0. Then the limit linit_>oo V'lO exists and is equal to 0. 

Proof. Consider the dynamical system 

c^Q 

in the phase space C x ]R. Choose ^ > such that the set W[lj^{l, cq) = {(1, c) | |c — co| < (5} is a 
local center manifold of the equilibrium (l,Co). 

By Theorem X.2.1 in [TT], we can choose e < S so small that every invariant set in an e- 
neighborhood of (l,co) is contained in the local center manifold of (l,co). ("00, c) has a backward 
orbit that is contained in an e-neighborhood of (l,c). Since ^{t) is and bounded with bounded 
derivative for all negative t, this backward orbit is pre-compact and so has a well-defined a-limit 
set. This a-limit set is invariant and contained in an e-neighborhood of (1, cq), and hence is a subset 
of Wfg^{l,co). Since c is constant, it follows that a(V'o,c) = (l,c). In particular it follows that 
limt_>._oo i^it) = 1, as desired. This completes the proof of the first statement. 

The proof of the second statement is similar. Choose 6 > such that the set Wf^^{Q,co) = 
{(0, c) I |c — Col < (5} is a local center manifold of the equilibrium (0,co). Since the w-limit set of 
(■00, c) is invariant and uniformly close to 0, it is a subset of W[^^^{0, Cq). Since c is constant, it follows 
that limt_j.oo ^{t) = 0, as desired. □ 



To complete the proof of part 4 of Theorem 1.4 it remains to show that 0c is in fact initially 



decreasing — that is, we have to rule out that 0c is not the "positive" branch of W^{1). We only 
sketch the argument; it is similar to that at the close of the proof of Proposition |4.1| If 0c is not 
initially decreasing, it is initially greater than 1. Given any 77 > 0, we can choose a;o such that 
< 0(x) — 1 < 77 for all X < xq and such that 

0c(a;o) > 0c(a;o - t), re [0, r]. 

Since g is on a neighborhood of 1, our hypotheses on the gradient of g at 1 imply that, if 77 is 
small enough, then 0c(a;o) > 0. 0'(x) will similarly be strictly positive, for x > xq, at least until 
0(x) = rj. In particular, 0c is not close to 0C(,, a contradiction. 



This completes the proof that U is open and hence completes the proof of Theorem 1.4 

5 Behavior at G dU 

We now turn to the proof of our second main theorem. 

Proposition 5.1. Assume that g is and that (G3) holds. For each c* G dU with c* > h{SIg{Q)) 
there are two globally hounded non-constant solutions to ip^ and -0^, neither of which are con- 
nections from 1 to 0, such that ip^{oo) — and ip^(—oo) ~ 1. 

Proof. Let € dU, choose {cfc} C U with c^ — !• and let 0^ denote the corresponding initially 
decreasing fronts with wave speed c^. The 0fc are uniformly bounded by some M by Lemma |l.3| By 
Proposition 2.5 a uniform-on-compact-sets limit 0* is a solution of ^ that is defined for all x e M. 



Since U is open and c* G dU by hypothesis, 0* does not connect 1 to — that is, 0* does not satisfy 
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both boundary conditions oo) = 1 and 0*(oo) = 0. We claim, though, that by translating the 
(pk properly we can arrange for <j)^ to satisfy either of these conditions. 

Since all of the (pk are initially monotone and decreasing until they cross 0, we can shift all of the 
(f)k to satisfy (^^(0) = 1/2, with (fik decreasing for a; < 0. Then satisfies these same conditions, 
and cxd) = 1. In this case, is the ip^ of the theorem. 



Recall that we are assuming that c* > 6(V.g(0)). By Proposition 2.7 then, there is some e such 
that for all large fc, the basin of attraction of for (and c*) contains the ball in C of radius e. 
Each <j)k approaches and so there is, for each fc, a final time when |(/)/j(a;)| = |. Translate each 
(pk to make this final time at time 0. Then (/)*(0) = ±e/2, and < | for all x e [0, 1]. By 

Proposition |2.7[ 0* must approach as a; — ?■ oo. 0* is not, however, identically zero. In this case, 

is the -0" of the theorem. □ 

6 Examples 

In section [l.l I we mentioned [27], in which Peletier and Rodriguez study the equation 

U'n = -Un + h{Un-l), h{u) = | ^ ^ j^^ (30) 

(we have changed variables to put (30) into the framework of equation ([2])). In our notation, for 
this equation c(V(;(0)) = c(/3) — c(— 1,2) « 4.31107. Using a mixture of numerical and analytical 
techniques, in [27] the authors explore the behavior of fronts as c decreases. In particular, they find 
values < c^nb < cuf < c(— 1, 2) such that initially decreasing solutions of Q with (j){—oo) — 1 

• approach monotonically for c > c(— 1, 2); 

• approach nonnionotonically for c G (c{,i/, c(— 1, 2)); 

• are bounded but do not approach for c G [c,j„f,, Cbif]] 

• are unbounded for c < Cunb- 

The value Cbif is where the leading roots of the characteristic equation cross the imaginary axis. For 
this example, therefore, in our notation we have — c(Vg(0)) and Cf = b{\/g{0)) — otherwise 
put, the bifurcation values Cm and c/ are exactly what one would expect given only the linearization 
of (|4| at 0. 

We mention again the work of Hsu, Lin, and Shen in [15', which includes a thorough discussion 
of the behavior of traveling wave solutions of various speeds for a particular family of lattice dif- 
ferential equations with unidirectional coupling. In the equations considered in |15] there are three 
equilibria fi^ < ijP < n^. For wave speeds large enough there is a monotone traveling wave solution 
connecting /i" to Under appropriate conditions, as the wave speed decreases the sequence of 



wave behaviors is reminiscent of that for ( 30 1 : we have monotone waves connecting fjP to /i"*" , then 
bounded nonmonotone waves oscillating about fjP a,s x ^ —oo, then (for a particular wave speed) 
a monotone wave connecting /i~ to /x"*", and then waves satisfying pioo) = /i+ that are unbounded 
below. Somewhat more generally, there is a range of wave speeds such that there is no connection 
from f/* to but such that waves are bounded. (See Theorem A in [15 ). 



Similarly, to return to the discussion of (30 1, the range c G [cunb, Cbif] corresponds to bounded 



traveling waves that do not connect 1 to 0. What is apparently happening here (expressed in terms 
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of our notation) is that, as c drops below Cf, the initiaUy decreasing branch of the unstable 



manifold at 1 (what is called, in the c = c/ case, in Theorem 1.6) is connecting 1 to a periodic 



orbit about 0. (Indeed, for (30), Peletier and Rodriguez find a subinterval of [cunb,Cbif] for which 
waves are eventually periodic, and there is an analogous interval of wave speeds for some of the 
class studied in [TJ]. Waves that are literally eventually periodic, though, arise as artifacts of the 
piecewise linear feedback functions used in both _27 and [15 .) It seems, though, that there should 
be examples of ^ for which initially decreasing solutions of (4j5) are unbounded for all c < cy, and 



examples for which ip'^ connects 1 to a nonzero equilibrium solution. 

Our goal in this section, broadly speaking, is to suggest the breadth of possible behaviors of 
solutions of (4|5) as c varies. We will confine our attention to the following specialized single-delay 



problem: 

c0'(a;) = (/.(x)-/(0(x-l)). (31) 

In this particular context, the assumptions (Gl) may be expressed as follows (recall section [lT| ) : / 
is a function such that 

. /(O) = 0, /(I) = 1, and f{s) ^ s for s e (0, 1); 

. /'(O) > 1 > /'(I) > 0; 

• There is some /3i > 1 such that s < f{s) < /3is for s G (0, 1). 
Note that, in this context, since g{so, si) = — sq + /(si) we have 

c(V5(0)) - c(-l, /'(O)) and 5(V.g(0)) - b{-l, /'(O)). 



Note also that, for (31), it is sufficient for (G3) to hold that / be bounded. 

We present the following three examples. The first two show that the linearization at does not, 
in general, provide exact information about the locations of Cm and c/. 

A: Given /'(O) and z > 0, / can be chosen so that c„j > c(— 1, /'(O)) + z. 

B: / can be chosen so that c/ lies strictly between 6(— 1, /'(O)) and Cm- 

C: Non- initially decreasing fronts can indeed exist. 

6.1 Example A 

Let a > 1 and z > be given and let 7 be so large so that c(— 1,7) > c(— l,a) + z. We will define 
a two parameter family of functions fr.rj which satisfy (Gl) as follows. Choose and fix functions 
/o.o and /i.i which satisfy (Gl) with /oo(0) = a and /{ ^(0) — 7. Moreover, let /o,o and /i^i be 
concave down so that a and 7 also play the role of /3i for /o,o and respectively. It follows from 



Theorem 1.4 that for c > c(— l,a) the lattice equation ([31)) with / given by /o,o has a monotone 



front connecting 1 to 0. On the other hand, for c < c(— 1, 7), the equation (31 ) with / given by 



has no monotone front connecting 1 to (Proposition 3.5). 

For (t, ?/) e (0,1] X [0,1], let pr.ri be a smooth cutoff function which takes values in [0,1], is 
increasing, and satisfies 



Pr.rjis) 



s<t/2; 

77 S > T. 
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Define a two-parameter family of feedback functions fr.ri by 

fr,r,{s) ■■= (1 - Pr,^(s))/o,o(s) + Pr.r, (s)/l,l («) 

— that is, f-r,r] is equal to /o,o for s < r/2 and equal to the convex combination "rj of the way from 
/o,o to for s > T. Such an fr^r] is illustrated in Figure 1. 



Figure 1 




s 



Note that /^ ,,(0) = a for (r, 77) G (0, 1] x [0, 1]. When ?/ = we have fr,o = /o,o and thus (jSTj) 
with / = fr.o admits a monotone front for c > c(— 1, a). 

We claim the following: For each c e [c(— 1, a), c(— 1, 7)) — in particular, for some c > c(— 1, a)+2; 



there is some and some 77 e [0, 1] such that (31 1 with wave speed c and / = /r,,r) admits no 



monotone front connecting 1 to 0. For otherwise, an argument similar to that in the proof of 
Propositions 2.5 and 3.4 would imply that (31 1 with wave speed c and / — fi^i admits a monotone 



front connecting 1 to 0; this contradicts Proposition 3.5 (The proof of Proposition 2.5 has to be 
modified to apply to a convergent sequence of feedback functions rather than a convergent sequence 
of wave speeds.) 

The point is that /'(O) alone provides only a lower bound on c„j; in general, all we know about 
Cm is that it hes somewhere in the interval [c(— 1, /'(O)), c(— 1, If / is concave down, of course, 
since /3i can be taken equal to /'(O) we necessarily have = c(— 1, /'(O)). 

6.2 Example B 

In this example we show that c/ can be strictly greater than 6(Vg(0)) = &(— 1,/'(0)). Similarly to 



the last example, the point is that the linearization of (31 1 at does not determine the range of wave 



26 



speeds for which fronts exist. (We remark that, while it seems plausible both that U is connected and 
that U has lower bound 6(— 1, /'(O)), we have not proven either of these claims.) Since our examples 
will use smooth bounded feedback functions /, hypothesis (G3) holds and c/ > 6(— 1,/'(0)) implies 
that the solutions '4>'^ and ■0° described in Theorem 1.6 exist. 



Wc first consider equation (311 with feedback function f = fe, where e £ (0,1/2) and /c is a 



smooth, odd, increasing function satisfying 

Ms) = 



2s, s€[0,6] 
1, s > 2e. 



For definiteness, let us take c = 5. Observe that that 5 > c(-l,/'(0)) = c(-l,2) > &(-l,2). We 
show that, for e small enough, f — fe, and c ~ 5, the initially decreasing branch 0_ of the unstable 



manifold of (31 1 at 1 coincides, after finite time, with a nontrivial periodic solution oscillating about 
zero. 
Define 

S {V^ G C : I : ip{0) = 2e and iP{x) > 2e for x £ [-1, 0]}, 

where the notation is chosen to be suggestive of a Poincare section. Define E' C E to be the subset 
of '0 S S such that V't G ^ for some t > 1. Define the first return time t : E' — )■ [l,oo) by 
r('0) = mit>i{t : I -04 e E}. Define the first return map P : E' E' by P(V') = V't(i/>)- Fixed points 



of P correspond to periodic orbits of (31). We claim that P has a fixed point 0*; moreover, the 
initially decreasing branch of the unstable manifold at one, which we now denote by 0, coincides, 
after a finite time, with 0*. The two observations, (i) P takes at most one value and (ii) € E', 
imply the existence of a fixed point 0* = P(0) and the fact that 0( is eventually periodic. To see 
that P cannot take more than one value, suppose that 0"^ and -0^ are both in the range of P. Let 



'0^(i) and 0'^(i) denote their forward continuations under (31) with f = ft- Since is constant for 
x>2e and tp^,ip'^ £ E, it follows that -0^ and -0^ agree on [0, 1]. Since ?A^(t) and ip'^{t) agree on an 
interval of length one, it follows that they agree for all t > 0. 

We now establish that 0_ G E'. Observe that (a particular translate of) is given by 

0(x) = (2e-l)e^/5 + l, x<0. 

Note that 0o € E. To show that 0o £ E', we claim that there is some > such that 0_ (x*) = — 2e 



and 0-(a;) < — 2e for all x £ [x* — 1, x*]; since the feedback in (31 ) is odd when / = /<;, by symmetry 
we have that 0_(2x*) = 2e and 0-(x) > 2e for all x £ [2x* — l,2x*]. Thus the translation of 
02x. £ S and hence 0o £ E'. 

We now prove the claim. For x £ [0, 1] 0_ satisfies the ODE 

f X (t>~ix) - 1 
= ^ 

and so is still given by 0_(x) = (2e - l)e^/5 + 1. Therefore we have, for e small enough, that 
0_(x) — — 2e for some x — xq £ (0, 1) — in particular, 

0_(x) = -2e for Xq := 5 In Q ^ < 1. 

We also have 0-(l) = (2e — l)e^/^ + 1 < — 2e. For e small, xq is small and 0-(l) is approximately 
—0.22. We wish to show that —1 < 0-(x) < — 2e for all x £ [xq,xo + 1], for then 0'_(xo + 1) = 
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{4>^{x() + 1) + l)/5 > 0, and 0_ will continue to increase until the time in our claim. For e small, 
4>-{x) is certainly in the desired range for x £ [xq, 1]. Since /<= is increasing, for x G [1, 1 + xq] we 
have that — 1)) is between —1 and 1. Thus, crudely, for x e [1, 1 + xq] we have 

5 ^^-(-)^ 5 ■ 

By taking e small enough and comparing (x) to the sub- and supersolutions suggested by the 
above inequality on the small interval [1, 1 + xq], we obtain the desired result. Thus, even though 
5 > &(— 1, /^(O)), 5 is clearly less than c/ when f = fe and e is small enough. 

We remark briefly that we have used a very specialized that guarantees that P takes only one 
value in order to simplify our analysis. However, we believe the phenomenon of fronts limiting to 
connections between equilibria and periodics to be more robust. An argument similar to that used 
above can, in many cases, show that e E' and hence that the initially decreasing branch of the 
unstable manifold exhibits recurrent dynamics. The presence of a periodic solution as a potential 
w-limit set for cj) can, in many cases, be shown to exist via Schauder's Theorem together with the 
fact that P is compact. 

Wc now turn an example of a different kind of behavior as c reaches and drops below c / (we have 



not rigorously proven what we are about to describe). Consider equation (31 1 with / = smooth, 
increasing, and satisfying the following, where a > and e G (0,a/2): 

r 1, s> 1/2; 

Us)^\ 2s, ha/2 + e, 1/2]; 

[ -a, s< -a/2. 

We conjecture that a and fa can be chosen so that the following hold: 

• c/ is as close to c„i = c(— 1, 2) as desired; 

• At c — Cf, 4'{x) — —a for all sufficiently large x; 

• For c < Cf, 0(x) — — oo as X — )■ oo. 

In this example, therefore, connects two equilibria (and, borrowing notation from f27i, Cf,i/ = 



There is clearly an equilibrium solution at —a. Any solution of (31 1 with 4>{xq) < ~a and 
4>{x) > 0(xo) for all x < xq will satisfy (/"(oo) = —oo. Also, as c — J- c/ from above, the minimum 
value of the front <^c connecting 1 to will decrease continuously in c. If c/ is small enough to 
allow this minimum value to approach —a, then at c = Cf this minimum value is equal to a and 
(pcfix) = —a for all sufficiently large x, and c = c/. 

We have not proven rigorously that c/ is actually low enough for this to work — that is, we 
have not ruled out a narrow parameter range for which 0c (x) > —a for all x but approaches a 
periodic solution rather than the equilibrium at 0. Numerical evidence with a = 1/4 and e = 0, 
though, suggests that c/ « 2.1 and that for all c < c/ the initially decreasing branch of the unstable 
manifold at 1 is unbounded. Note that 2.1 is well above b{—l, 2) « 1.65. We could, in fact, make c/ 
as close as we like to c„i — c(— 1, 2) by taking a close to 0. 

Figure 2 shows the initially decreasing branch of the unstable manifold at 1 for various values of 
c near 2.1. Note that, on compact subsets, the solutions that ultimately approach appear to be 
converging to a solution that approaches —a — that is, the solution ■0^. 
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Figure 2 
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6.3 Example C 

While our focus in this paper is on initially decreasing fronts, we emphasize that non-initially de- 
creasing fronts can exist: in the notation of Lemma 2.4, (f)^(oo) = is possible. We present a 
concrete example, somewhat informally. 



We consider (31) where / is an increasing function satisfying 



m = 




s < 1/2 -e; 
se [1/2,3/2]; 



> 3/2- 



where e > and such that we can take /3i slightly larger than 2. We take c — 5 > c(— 1, /3i)); in this 
case there is a decreasing front Since /(s) = 1 for s near 1, for large negative x any solution of 



(31) with 0(— oo) — 1 satisfies the ordinary differential equation 



and is therefore of the form 



<P'{x) 



(f>{x) = Ce^/^ + 1. 



Taking C < yields the initially decreasing branch ^_ of the unstable manifold, and taking C > 
yields the initially increasing branch of the unstable manifold. In particular let us take the 
function (f>^ defined by 



(j>+{x) 



1 



,x/5 



1 



29 



for all X < 0. Observe that 0+(O) ~ 3/2. For a: > let us continue as a solution of the delay 



equation (31). We wish to show that (/)+(oo) = 0. We will do this by proving the following claim: 
there is some a;* > such that < 1 and (/)+(a;) G [1/2,3/2] for all 2; S [x* — l,a;*]. It follows 

that, for all a: > a;*, coincides with a translate of and so is a front. 

will continue to be of the form (l/2)e^/5 + 1 for x e [0, 1]. Observe that 0(1) = (l/2)ei/^ + 
1 « 1.61. Let us take e < 0(1) — 3/2 so that x\ = infja: > : 4>(x) — 3/2 + e } is less than 
1. Observe that, for all x G [l,a;i + 1], the quantity (\){x — 1) lies between 3/2 and 3/2 + e and so 
f{(f>{x — 1)) lies between 1 and 5. If e is small, we have that Xi is close to and that 0(1 + Xi) is 
close to 0(1). In particular, we may assume that 0(a;) > 3/2 + e for all x G {xi,xi + 1]. 
Let us write 

X2 = inf{ X > xi + 1 : 0(x) = 3/2 + e }• 
On the interval [xi + 1, 0:2 + 1], (pix) will satisfy the ODE 

[X) = — 1. 

In particular, 0(x2 + 1) will be equal to 

(0(a;2) - 5)6^/5 + 5 w -^e^/^ + 5 w 0.725. 

For X G [x2,a;2 + 1], we have (crudely) 0'(x) < —3; and for x £ [x2 + 1,X2 + 2], we have 0+(a;) > 
(0+(a;) — 5)/5. Therefore, in particular, for e small enough, 0+(a;) reaches 3/2 shortly after X2 and 
does not reach 1/2 too soon after X2 + 1- Thus there is some a;* > 3:2 + 1 such that 0+(a;*) < 1 and 
0+(a;) G [1/2,3/2] for all a; G [a;* — l,a;*]. This proves our claim. 

Figure 3 illustrates the two traveling wave solutions of the equation just described, with a small 

e. 
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Figure 3 
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